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ABSTRACT 

A two-dimensional particle-in-cell simulation is performed to investigate weakly magnetized perpen- 
dicular shocks with a magnetization parameter of a — 6 x 10 -5 , which is equivalent to a high Alfven 
Mach number Ma of ~ 130. It is shown that current filaments form in the foot region of the shock 
due to the ion-beam- Weibel instability (or the ion filamentation instability) and that they generate a 
strong magnetic field there. In the downstream region, these current filaments also generate a tangled 
magnetic field that is typically 15 times stronger than the upstream magnetic field. The thermal 
energies of electrons and ions in the downstream region are not in equipartition and their temperature 
ratio is T c /T; ~ 0.3 — 0.4. Efficient electron acceleration was not observed in our simulation, although 
a fraction of the ions are accelerated slightly on reflection at the shock. The simulation results agree 
very well with the Rankine-Hugoniot relations. It is also shown that electrons and ions are heated in 
the foot region by the Buneman instability (for electrons) and the ion-acoustic instability (for both 
electrons and ions). However, the growth rate of the Buneman instability is significantly reduced 
due to the relatively high temperature of the reflected ions. For the same reason, ion-ion streaming 
instability does not grow in the foot region. 

Subject headings: instabilities — magnetic fields — plasmas — shock waves — supernova remnants 
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1. INTRODUCTION 

A large volume of the universe (including interstellar and intergalactic space) is filled with hot, tenuous plasmas. 
Coulomb collisions between charged particles rarely occur in these plasmas and the plasma dynamics are dominated 
by collective phenomena involving particles and electromagnetic fields (e.g., plasma oscillations). Hence, such plasmas 
are known as collisionless plasmas. Even in collisionless plasmas, some kinds of "shocks" occur. These shocks generally 
have very complex formation mechanisms that involve various kinetic processes, including electrostatic instabilities, 
electromagnetic instabilities, and compression of background magnetic fields. The shocks driven in supernova remnants 
(SNRs) are considered to be collisionless ones. 

SNR shocks propagate in the interstellar medium, which has weak magnetic fields of typically ~ 3/xG. In the context 
of shocks in magnetized plasmas, the strength of the magnetic field, -Bo, is frequently expressed in terms of the 
magnetization parameter or the sigma parameter, which is defined as the ratio of the magnetic energy density to the 
bulk kinetic energy density of the upstream plasma (both are measured in the shock rest frame) . For nonrelativistic 
cases, it is given by 



n e0 (m c + mi)V s l/2 



= Ma 4 , (1) 



where n o is the electron number density in the upstream plasma, m c is the electron mass, mi is the ion mass, V^h is 
the shock speed, and Ma is the Alfven Mach number. For example, for the shock in SN1006 (except the North- West 
region), n C Q ~ 0.05 cm -3 and V s h ~ 4900 km s _1 were inferred ([Acero et al. Il200"7l ) so that <r~4x 10~ 5 . The recently 
disco vered 'youngest' SNR Gl.9+0.3 is considered to have a shock velocity of V s h ~ 14, 000 km s -1 (|Revnolds et al. I 
2008); assuming that n C Q ~ 0.1 cm -3 gives a ~ 2 x 10~ 6 . The sigma generally lies in the range 10~ 6 < a < 10 -3 
for shocks in young SNRs; these shocks are thus very low-cr shocks (or, equivalently, very high Alfven Mach number 
shocks). 

Magnetized shocks have been extensively investigated, especially perpendicular shocks in which the background 
magnetic field is perpendicular to the shock normal. The structure of perpendicular shocks in the supercritical regime 
(Ma > M^, where M A ~ 3) is known to some extent: a fraction of the incoming ions are reflected at the shock front 
(called the 'ramp') and the reflected ions form a slightly dense region, referred to as the 'foot', in front of the ramp. 
The ions also accumulate immediately behind the ramp and generate a strong magnetic field there, which is called the 
(magnetic) 'overshoot'. Over the last decade, one-dimensional (ID) particle-in-cell (PIC) simulations that can model 
the kinetic d ynamics of both electrons and ions have been performed to investiga te high-Mach-number perpendicular 
shocks (e.g.. IShimada fe Hoshino 1120001 : ISchmitz et al. 112002^ iScholer et al. 1120031). Recently, several two-dimensiona l 
(2D) simulations have also been performed (e.g.. lUmeda et aL 

KQ8; A mano fe Hoshino Il200l iLembege et al. Il2009h . 

However, most simulations have been conducted for relatively strong background magnetic fields (a > 5 x 10~ 3 or 
Ma < 15). It is thus desirable to perform simulations for weaker background fields. 

On the other hand, it was recently demonstrated that certain kinds of collisionless shocks can occur even in un- 
magnetized plasmas at relativistic shock sp eeds by performing two- or three-dimensional (3D) PIC simulations (|Kato I 
l2007t ISpitkovskv 1120081 : IChang et al~ll2008D . In these shocks, the beam-Weibel instability (or filamentation instability) 
is driven in the transition region of the shocks between the counterstreaming electron-positron beams in pair plasmas 
or between the counterstreaming ion beams in electron-ion plasmas and generates strong magnetic fields there. These 
generated fields provide an effective dissipation mechanism for collisionless shock formation and are hence often referred 
to as ''Weib el-mediated shocks." [Note that the beam-Weibel instability is driven by the counterstreaming be ams (c.f. 
Fried IH959I ) an d it differs from t he ord inary Weibel instability, which is driven by a temperature anisotropy ()Weibel I 



19591 ) fsee also iDavidson et al. I (fl972l)).1 These shock s can be driven b y relativistic phenomena, such as gamma-ray 



bursts and their afterglows jlMedvedev fe Loeb I Il999t iBrainerdl 120001 ) jets from active gala ctic nuclei, and pulsar 



winds (jKazimura et al. Ifl998h . It was also shown in our previous paper (jKato fe Takabe l l2008) that this kind of shock 
can form in unmagnetized electron-ion plasmas even at nonrelativistic speeds. The beam-Weibel instability can also 
be important in weakly magnetized shocks. As was shown in our previous paper, the magnetic field generated by the 
ion beam-Weibel instability reaches a few percent of the upstream bulk kinetic energy and this value is much higher 
than the background magnetic field around typical SNR shocks. Therefore, the ion beam-Weibel instability may play 
an important role in the formation of weakly magnetized shocks. Magnetized shocks have been extensively investigated 
by ID simulations. However, ID simulations cannot consider the beam-Weibel instability because its wave vector is 
perpendicular to the flow direction. Therefore, it is essential to perform multidimensional simulations to investigate 
the formation process of weakly magnetized nonrelativistic shocks. 

Collisionless shocks are also considered sites of particle acceleration. In particular, cosmic-rays with energies below 
10 15 eV are considered to be accelerated in SNR shocks. Inde ed, recent X-ray observations revealed that electrons 
are a ccelerated to energies of ~ 10 14 eV in several young SNRs (jKovama et al. I IT995: Lon g et al~ll2003l : iBamba et al. I 
2003) . It is widely accepted that first-order Fermi ac celeration or diffusive shock acceleration is the acceleration mech- 
anism (e.g.. lDru rv 1983; Blandf ord fe Eichler 1119871 ). However, it is currently not possible to determine the fraction of 
thermal plasma particles that are injected into the diffusive shock acceleration process (this is known as the injection 
problem). For electron injection in quasi-perpendicular shocks, the shock surfing accelera t ion has been investigated 
as an injection mechanism or even as an efficient acceleration mechanism (Hoshino 2001; McCl ements et al~l 120011 : 
Hoshino & Shimada 2002) . Howe ver, several researchers have recently shown that the shock surfing accelera tion pro- 
cess is in fact inefficient in 2D fe.g.. lDieckmann fc Shuklall2006t lOhira fc Takaharall2007l : lUmeda et al. 1120081) and that 
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it is efficient in ID because of the symmetry of the system. Instead, lAmano fc Hoshino I (|2009T > showed that another 
acceleration process can operate in 2D in which a fraction of electrons are reflected in the foot region by small-scale 
electrostatic waves generated by the Buneman instability. They are then accelerated by the motional electric field as 
well as being directly accelerated by the electric field when they resonate with the electrostatic waves. Electrons can 
be accelerated up to about the upstream ion bulk energy by this mechanism. Thus, multidimensional effects can play 
an essential role in the acceleration/injection mechanism. 

In addition, collisionless shocks can be sites of magnetic field amplification/generation. Recent X-ray observations 
suggest that ma gnetic fields of the order of hundreds of microgauss or even milligauss may be generated in the vicinity 
of SNR shocks (|Vink fc Lam ing 20031 FVolk et al. 1120051 : IUchivama~l l2007t ). Several mechanisms have been proposed 
for this m agnetic fiel d amplification, including a nonresonant instability dri ven by high-energy particles accelerated 
in sh ocks (|Bell ll2003 ) and magnetohydrodynamic turbulence behind shocks (jGiacalone fc Jokipii 1 120(371 : llnoue et al. I 
2009). The mechanism may be related with the microscopic kinetic processes associated with shock formation itself; 
it should in principle be possible to investigate this by performing large-scale PIC simulations. 

In this study, we investigate the formation and structure of perpendicular shocks for very low a and the pro- 
cesses responsible for particle acceleration and magnetic field generation by performing 2D PIC simulation, which 
can appropriately model the beam-Weibel instability. Because of the capability of the computer, we used a reduced 
ion-to-electron mass ratio and a shock speed (V s h ~ 0.3c, where c is the speed of light) that is much higher than 
realistic ones for SNRs (V s h ~ 0.01c) in the simulation. 

2. METHOD 

We investigated collisionless shocks in electron-ion plasmas with weak background magnetic fields by performing a 
2D PIC simulation. The simulation code is a relativistic, electrom agnetic, PIC code with two spatial and three velocity 
dimensions developed based on a standard method described by iBirdsall fc Langdon I (|1991h . The basic equations of 
the simulation are Maxwell's equations and the (relativistic) equation of motion for particles. In the following, we 
regard the simulation plane as the x — y plane and we take the z-axis to be perpendicular to the plane. We take r = lu^ 
to be the unit of time and the electron skin depth A c = cujp^ to be the unit of length, where u; pe = (47m oe 2 /m ) 1 / 2 is 
the electron plasma frequency defined for the electron number density in the far upstream region, n C Q. The units for 
electric and magnetic fields are E* = £>„ = c(47rn e om e ) 1 / 2 . 

In the simulation, a collisionless shock is driven according to the "injection method." There are two walls, one on the 
left-hand side (smaller x) and the other on the right-hand side (larger x) of the simulation box and these walls reflect 
particles specularly. Initially, both electrons and ions are loaded uniformly in the region between the two walls with a 
bulk velocity of V in the +x-direction. The electrons and ions have equal temperatures in the upstream region. In the 
early stages of the simulation, particles near the right wall are reflected by the wall and then interact with incoming 
particles (i.e., the upstream plasma). This interaction generates some instability and eventually a collisionless shock 
forms. The frame of the simulation is the rest frame of the shock downstream; the shock propagates from right to left 
in the downstream rest frame. 

We consider a perpendicular shock in this paper; the initial magnetic field, Bq, is in the y-direction (i.e., in the 
plane) and its strength is determined by the sigma. However, since the shock speed is unknown before performing the 
simulation, in the following, the sigma is defined in the simulation frame with an upstream bulk velocity V instead of 
the shock speed V^h as 

- _ gg/g7[ , 

a = 7 1 u?2 /o ' ( 2 > 

however, the difference between these two sigmas is not large. With this definition of the sigma, the magnetic 

1/2 

field strength in the simulation frame is given by B$ = [(1 + nii/m c )d] (V/c)B*. The initial electric field, Eq, is 
determined so that it vanishes in the plasma rest frame (i.e., the upstream frame); this requirement causes the motional 
electric field in the simulation frame, Eq = —VBq/c, in the z-direction. The boundary conditions for both the particles 
and the electromagnetic field are periodic in the y-direction. 

3. RESULTS AND ANALYSIS 

We performed a simulation for a sigma of a — 10~ 4 . As mentioned above, we use a reduced ion mass of m\ — 30m o 
and a bulk velocity of V = 0.25c. The grid size is N x x N y = 16384 x 1024 and there are ~ 40 particles per cell per 
species. The physical dimensions of the simulation box are L x x L y = 3200A o x 200A o and thus the size of a cell is 
A.t = Ay ~ 0.2A C . The electron and ion temperatures are equal and are given by k-^T ] /m e c 2 = 1.25 x 10~ 3 , where 
fee is the Boltzmann constant. The thermal velocities are thus given by a = (2fcB? 1 /m c ) 1 / 2 = 0.05c for the electrons 
and ai = 9.13 x 10 _3 c for the ions. For these parameters, we have lu cc /lu pc = Bq/B* ~ 1.4 x 10~ 2 , the Alfven speed 

va ~ 2.5 x 10~ 3 c (thus, Ma = V/va — 100), and the plasma beta (3 ~ 26 (i.e., it is a high-beta plasma). The Larmor 
radii of the electrons and ions calculated for the background field and the upstream bulk velocity are r 9ie = 18A e and 
fg.i = 535A e , respectively. 

3.1. Overall structure 

Figure Q] shows the time evolution of the ion number density averaged over the y-direction. The shock transition 
region, or the "shock front" , appears as a steep increase in the number density. The shock structure and its propagation 
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speed abruptly change around ui pe t ~ 3000. This is because the shock structure undergoes a transition from an 
unmagnetized shock to a magnetized shock. Indee d, the structure for u u e t < 2000 is essentially the same as those of 
Weibel-mediated shocks in unmagnetized plasmas ([Kato &: Tak abc 200g), as discussed below. The transition time is 
of the order of the gyration time of the ions in the background field, T g = 27r/a; c i, where oj c ; = eBo/rriiC is the ion 
cyclotron frequency. For the unmagnetized shock (1500 < 0J pe t < 2000), the shock speed measured in the downstream 
frame is V^h.d ~ —0.16c. For the magnetized shock (uj pc t > 4000), it becomes I4h,d ~ —0.08c and that in the upstream 
frame and the Alfven Mach number are given by V^ n ~ —0.33c and Ma ~ 130, respectively. Thus, the sigma defined 
for the shock velocity is given by a = 5.9 x 10 -5 in this case. The shock speeds obtained here may have small 
uncertainties because they were obtained by eye-fitting the figure and also they may not be in the steady state yet. 
Since the formation of an unmagnetized shock is a consequence of the initial conditions and we are interested in the 
magnetized shock in this study, we mainly focus on the magnetized shock below. We discuss the unmagnetized shock 
at the end of this section. 




e 



Fig. 1. — Time development of the ion number density averaged over the ^-direction for the simulation with an upstream bulk velocity 
of V = 0.25c and a sigma parameter of a = 10 — 4 . The color indicates the number density normalized by that in the far upstream n e Q. The 
horizontal and vertical axes respectively represent x (in units of electron skin depth A e ) and time (in units of the electron plasma time Upe). 
The shock structure and its propagation speed change abruptly around uj pc t ~ 3000 due to the transition from an unmagnetized shock to a 
magnetized shock. The arrows indicate the obtained shock speeds for the unmagnetized shock (V^h d ~ —0.16c) and the magnetized shock 
(Vsh.d ~ -0.08c). 

Figure Q] shows that the shock wave almost reaches steady state near the end of the simulation (uj pc t ~ 8000). Figure 
[2]shows the ion number density at uj pc t — 8000. (Hereafter, we discuss the results at this time unless otherwise stated.) 
The upstream plasma flows from left to right and moves through the transition region (2350 < x/X c < 2550) and then 
reaches the downstream state. (The structure in Fig. Q] in x > 2700A o is an artifact due to the boundary and so in 
the following we discuss the structure in x < 2700A o .) Note that there are filamentary structures, which cannot be 
observed in ID simulations, in the upstream leading edge of the shock transition region (x ~ 2400 A c ). The filament 
radius is typically approximately equal to the ion inertial length, wh ich is the same as those in the "Weibel-mediated" 
shocks in unmagnetized electron-ion plasmas (jKato fc Takabe I l2008l ). Then, behind them, there is a highly fluctuating 
high-density region. In the downstream region (x > 2550A o ), the number density becomes almost homogeneous. 

Figures [21a) and (b) show profiles of the ion number density and the magnetic field strength averaged over the 
y-direction, respectively. It shows that the number density increases rapidly in the transition region and after the 
transition region it approaches ~ 4 times the upstream value. Figure [3](b) shows the root mean square of each 
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Fig. 2. — Ion number density at uj pe t = 8000. The horizontal and vertical axes represent x and y, respectively. The left- and right-hand 
sides are upstream and downstream of the shock, respectively. Many filamentary structures can be seen in the shock transition region 
(2350 < x/X c < 2550). 

magnetic field component together with the total magnetic held strength. This structure is similar to the well-known 
structure of supercritical perpendicular shocks in ID; the 'ramp' is at x ~ 2400A G and there is an extended 'foot' 
region in x < 2400 A e as well as an 'overshoot' region in 2400 < x/X c < 2470. It is evident that strong magnetic helds 
are generated in both the shock transition region (or the overshoot) and the downstream region. The energy density 
of the magnetic held reaches ~ 15% of the upstream bulk kinetic energy density (measured in the downstream rest 
frame) in the shock transition region and ~ 2% in the downstream region. It is notable that B x and B z , which are 
generated by the current filaments of the ion-beam- Weibel instability (see below), are comparable with B Vl which is 
mostly generated by the upstream background field. These B x and B z helds as well as the B y field contribute to the 
dissipation of the shock. Since these filaments and the magnetic field are never generated in ID simulations, the shock 
structure may differ significantly from those in ID cases. 



8-r 




Fig. 3. — Profiles of (a) the ion number density normalized by the upstream density n c o and (b) the root mean square of each magnetic 
field component (B x , blue curve; B y , red curve; B z , green curve) and that of the total strength dashed black curve), where all the 

components are normalized by the upstream background field Bo. 

Figure [4] shows phase-space plots of the electrons and the ions. Here, each component of the four velocities (uj = 
jVj/c, where j = x,y,z and 7 = (1 — w 2 /c 2 ) -1 / 2 is the Lorentz factor of the particle) are plotted as a function of 
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the x-coordinate. Both electrons and ions from upstream are mostly dissipated and isotropically thermalized through 
the transition region (2300A o < x < 2450A o ). It is observed that a fraction of ions are reflected at x ~ 2400A C (i.e., 
the ramp) and then gyrate back downstream with slight acceleration forming the foot structure . This is a well-known 
chara cteristic of supercritical shocks and has been observed in many numerical simulations fe.g.. lLerov Ill98ll iBurgess I 
1989). In contrast, the electrons have no prominent substructures in phase space. 




2200 2300 2400 250O 2600 27O0 2200 2300 2400 2500 2600 

x/ A x. / "k 



Fig. 4. — Phase-space plots of electrons (left panel) and ions (right panel). The (from top to bottom) x, y, and z components of the 
four velocities are shown in each panel. Both species are mostly thermalized within the shock transition region. A fraction of the incoming 
ions are reflected at x ~ 2400A e . 

3.2. Foot dynamics 

Figure [5] shows the distribution of the ion kinetic energy measured in the upstream frame, i^kin.u = (ju — l)mc 2 , where 
7 U is the particle Lorentz factor measured in the upstream frame, as a function of the x-coordinate (in the downstream 
frame) at uj pc t = 8000. In this figure, the incoming and reflected ions in the foot region (2200 < x/X c < 2400) can 
be clearly distinguished from each other using a threshold energy of, for example, E^^/nicC 2 = 1: the reflected ions 
with i?kin,u/™c c2 > 1 an d the incoming ions with E^i m ^ n /m e c 2 < 1. The reflected ions are further divided into two 
populations: those streaming upstream measured in the shock rest frame and those streaming downstream. Thus, we 
can investigate the foot dynamics on the basis of a simple fluid model that consists of a single electron fluid and three 
ion fluids (inc oming i o ns, re flected ions streaming upstream, and those streaming downstream), which is similar to the 
model used in ILerovH ()1983l ). For convenience, we denote the electrons, the incoming ions, the reflected ions streaming 
upstream, and those streaming downstream by the symbols e, I, R_, and R+, respectively. 




Fig. 5. — Kinetic energy of the ions measured in the upstream frame. 

Figures [Sfa) and (b) show the mean velocity of each fluid component in the x- and z-directions and Fig.|5Jc) shows 
the number densities normalized by the upstream number density. In the downstream region (x > 2400A C ), only the 
values for all the ions are shown because classifying the ions by the above method is meaningless in that region. It 
shows that the mean velocities and the number density of all electrons (thick curves) and for all ions (dashed white 
curves) agre e well with ea ch other in both the upstream and downstream regions, indicating that the massless electron 
fluid model (|Lerov II1983D holds well, at least on average, even in this high Mach number and low ion-to-electron mass 
ratio case. 
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3.2.1. Electrostatic instabilities and heating 

The local temperatures of the respective components were calculated using the mean velocities obtained above and 
they are plotted in Fig. [7] It shows that in the foot region, the electrons are heated in two steps: the first step in 
x < 2350 A c (region 1) and the second step in 2350 < x/X c < 2400 (region 2). The incoming ions are also heated in region 
2. This sequential electron heating process t ogeth er with the ion heating suggests that the model for very high Mach 
number shocks proposed by iPapadopoulosI (|1988f) is valid in the foot region, in which the incoming electrons are first 
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heated by the Buneman instability (Buneman 1 fl958h for reflected ions (|Auer et al. 1119711) and subsequently, after the 
electrons have been heated to a certain temperature, they are further heated by the ion -acoustic instability for incoming 
The latter instability can also heat the ions. This process has been studied by|Cargill & Papado poulos I (fl 988) 



ions. 



for M A ~ 
(200(3) for 



50 and 500 with hybrid simulations with a phenomenological resistivity and also by iShimada fc Hoshino I 
M A - 10.5 with ID PIC simulations. 
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Fig. 7. — Temperature profile of each component. Both incoming electrons (thick solid curve) and incoming ions (thin solid curve) are 
heated in the foot region (x < 2400A e ). 



This heating process is expected to operate in very high Mach number shocks and so it should also operate in the 
present case (Ma ~ 130). Figure [5] shows several quantities of each ion component obtained from Figs. [5] and [7J 
specifically, it shows profiles of the mean velocity relative to the electron velocity in the x-direction, the electron-to- 
ion temperature ratio, and the number density normalized by the local electron number density n e (x). In region 1 
(x < 2350A o ), the reflected ions streaming upstream have a significantly higher velocity relative to the electrons than 
the electron thermal velocity. On the other hand, in region 2 (x > 2350A e ), the electron-to-ion temperature ratio for 
the incoming ions increases to a large value and also the incoming ion velocity relative to the electron velocity becomes 
large due to the large deceleration of the electrons [see Fig.[5][a)] so that it becomes higher than the ion-acoustic speed, 
c s = (feeTc/mi) 1 / 2 . These conditions are indeed preferable to the Buneman instability in region 1 and the ion-acoustic 
instability in region 2. 

Here, we show the instabilities that operate in the foot region by performing local linear analysis with the fluid quan- 
tities (namely, the mean velocities, the number densities, and the temperatures shown in Figs. [5] and [7|). Approximating 
the distribution of each component as a Maxwellian distribution, 



f(> (Vx,Vy,V z ) = 



■ exp 



(v x - V StX ) 2 + vl + (v z - V s , z f 



(3) 



where s = (e, I, R_, R+), n s is the number density, V SiX and V s<z are the streaming velocities in the x- and z-directions 
respectively, and a s = (2kBT s /m s ) 1 ^ 2 is the thermal velocity. We solve the following dispersion relation for the 
electrostatic mode with the wavevector in the x-direction: 



fc 2 + 2^fc 2 s (l + C s ^(Cs)) =o, 



(4) 



where 



0J K 



( 47m s g s 2 



1/2 



(5) 
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Fig. 8. — Profiles of the quantities in the foot region for incoming ions (thin solid curves), reflected ions streaming upstream (dashed 
curves), and reflected ions streaming downstream (dot-dashed curves), (a) The velocities relative to the electron velocity in the ^-direction 
(absolute values). The electron thermal velocity (thick solid curve) and the ion-acoustic speed (thick dotted curve) are also shown, (b) 
The electron-to-ion temperature ratios, (c) The number densities normalized by the local electron number density, n e (x). 
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(6) 



The function Z(Q is the plasma dispersion function ([Fried fc Conte I [l 96 ID defined by 



OO 



Z(0 EE 7T 



-1/2 



e 



dz. 



(7) 



— OO 



TableQ]sumniarizes some quantities used in the following analysis. We use here the dispersion relation for unmagnetized 



TABLE 1 

Quantities obtained from the simulation. The units for the ^-coordinate, the velocity, the number density, and the 

TEMPERATURE ARE A , C, n o 0, AND m e C 2 , RESPECTIVELY. 



X 






n c 


To 






"i 


Tj 




V R -, Z 


n R _ 


Tr_ 






™R+ 


Tr + 


2300 


0.20 


-0.042 


1.1 


2.5 x 10~ 3 


0.24 


3.3 x 10~ 3 


0.99 


1.3 x 10~ 3 


-0.14 


-0.27 


0.069 


0.090 


4.6 x 10~ 3 


-0.39 


0.084 


0.14 


2350 


0.11 


-0.089 


1.6 


8.6 x 10~ 3 


0.23 


0.013 


1.0 


1.8 x 10~ 3 


-0.15 


-0.19 


0.37 


0.12 


5.5 x 10~ 3 


-0.34 


0.26 


0.19 


2375 


0.062 


-0.093 


2.1 


0.033 


0.22 


0.023 


1.0 


4.6 x 10~ 3 


-0.16 


-0.15 


0.68 


0.14 


0.013 


-0.32 


0.38 


0.22 



plasmas given by Eq. (Q} instead of that for magnetized plasmas because the magnetic field is sufficiently weak in the 
present case; the condition for the unmagnetized approximation is given by k 2 3> k 2 s where k cs = \/2\uj cs \ / a s ; in other 
words, for all species, the wavelength is much smaller than the Larmor radius defined for the thermal velocity. As 
shown below, the wavenumbers of the instabilities are typically kc/uj pc > 3 and k cs c/io pc is ~ 0.5 for electrons, ~ 0.15 
for incoming ions, and ~ 0.01 for reflected ions in the foot region. Therefore, the unmagnetized approximation can be 
used in this case. 

When performing the linear analysis with the local quantities at x = 2300 A c , we found an unstable electrostatic 
mode whose wavenumber (k x ~ 3.5w pe /c) and frequency (u/ ~ 0.048w pe in the rest frame of the reflected ions arc 
similar to those of the Buneman instability between electrons and reflected ions streaming upstream (k x ~ 2.9ui pe /c 
and u/ ~ 0.05wp C in the rest frame of However, the obtained growth rate (7 ~ 0.0057w pc ) is one order of 

magnitude smaller than the typical growth rate of the Buneman instability (7 ~ 0.087w pc ). This is because of the 
relatively high temperature of the reflected ions streaming upstream, Tr_ , as is shown in Fig. while the ordinary 
Buneman instability assumes that both species are cold. Figure [9] shows the maximum linear growth rates of this 
mode together with their wavenumbers calculated for the quantities at x = 2300A C while varying T^_. When Tr_ 
approaches zero, the growth rate becomes large and approaches to a typical value for the Buneman instability. Thus, 
we regard this mode as a Buneman instability between the electrons and the reflected ions streaming upstream with 
a reduction in the growth rate due to the relatively high temperature of the reflected ions. 

On the other hand, we found another unstable electrostatic mode at x ~ 2350A o . This mode has a maximum growth 
rate 7 max ~ 0.02w po and a frequency uj ~ 0.52cj pc at k x ~ 4.85w po /c. This leads to a phase velocity of ~ —0.017c in 
the in coming ion rest frame. For the same parameters, the dispersion relation of the ion-acoustic instability (Ichimaru 
1973) between electrons and incoming ions gives 7 max ~ 0.015cl> pc at k x ~ 4.0oj pc /c and a phase speed of c s ~ —0.016c 
in the incoming ions rest frame. Both agree well with each other and thus we regard this mode as an ion-acoustic 
instability between electrons and incoming ions. 

Figure flOl summarizes the results for this local linear analysis over the foot region. The Buneman instability develops 
upstream of the foot region (x < 2325 A e ), whereas the ion-acoustic instability dominates downstream of the foot region 
(x > 2325 A ). This feature is consistent with the evolution of the electron and incoming ion temperatures shown in 
Fig. [7] the electrons are first heated by the Buneman instability and then both electrons and incoming ions are heated 
by the ion-acoustic instability. Note that there is a region where both instabilities can coexist (2312 < x/X c < 2325). 

Since both the Buneman and the ion-acoustic instabilities are electrostatic modes, they are always associated with 
the charge density p and can be investigated through it. Figures EJa) and (b) show the charge density and its power 
spectrum in two rectangular areas in the foot region, namely x ~ 2300A o (where the Buneman instability dominates) 
and x ~ 2350A G (where the ion-acoustic instability dominates), respectively. The peak positions of these power spectra 
agree well with the wavenumbers for the maximum growth rates obtained by the linear theory shown in Fig. 1101 
Note that both spectra are not concentrated on the fc^-axis but extend in the fc y -direction. This results in the wavy 
appearance of both modes in real space (left panels) and is a well-known characteristic of both instabilities in multiple 
dimensions. 



As mentioned above, the ion number density in the foot region (Fig. [2]) contains many filamentary structures. 
Figure [12] shows that these filaments are associated with current filaments and filamentary magnetic fields. These 
filaments are similar to those observed in unmagnetized shocks, which are generated by the beam-Weibel instability. 
In our previous papers, we showed that the ion beam-Weibel in stability develops and generates currents filaments 
even for nonrelativistic flow speeds ()Kato fc Takabe 1 120081 12010D . Therefore, it is plausible that these filaments are 
generated by the ion beam-Weibel instability. To confirm this, we performed linear analysis in the same manner as that 



3.2.2. Filamentary structures 
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Fig. 9. — (a) Maximum linear growth rates of the electrostatic mode and (b) the wave numbers at the maximum growth rates as functions 
of the temperature of the reflected ions streaming upstream, Tr_, calculated for the quantities obtained from the simulation at x = 2300A C 
(except Tr_). Those for the original value of Tr— in the simulation are indicated by the arrows. 



used to obtain Fig. [10] except that we here consider the electromagnetic modes with wavevectors in the y-direction. 
In the present case, since the wavenumber is too low to employ the unmagnetized approximation for electrons, we 
solve the following dispersion relation in the electron rest frame, which includes the effect of the magnetic field for the 
electrons (the ions are assumed to be unmagnetized): 



where 
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Fig. 10. — (a) Maximum linear growth rates of the Buneman instability (filled circles) and the ion-acoustic instability (crosses) as 
functions of x. (b) The wave numbers at the maximum growth rates for the respective instabilities. 

Kz = ~\ (^) 2 ^° [Z{il) ~ Z{ ^- l)] + 2 ^ ( VsxVsz{l + asl (14) 



A. 



UJ 
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with 



ka e 

V sx = V SiX /a s , V e 



a s = t] s Z(t] s ), 



V S}Z /a s 



(15) 

(16) 
(17) 



In the above dispersion relation, the sums run only for the ion species, that is for s = I, R_ , R + . 

The results are shown in Fig. [13] by the solid curves. The mode is unstable in the foot region and it grows at a 
comparable growth rate to those of electrostatic modes (see Fig. [T0|) . The wavenumber obtained in the linear analysis 
near x = 2360A o is typically k y ~ 0.8; this value agrees well with the simulation result shown in FigfT2Td). Note that 
the real frequency of the mode (dotted curve) is zero; that is, it is a purely growing mode. 

This mode can be regarded as an ion beam-Weibel instability. Indeed, as shown in Fig. 1131 the maximum growth 
rate and the wave number essentially agree with those obtained from the dispersion relation for the beam-Weibel 
instability using the unmagnetized approximation: 



■2V S 2 (1 



0. 



(18) 



where V s is taken to be either V sx (shown by the dashed curves) or V sz (the dot-dashed curves) and s = e, I, R_,R+. 
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Fig. 11. — Charge density (left panels) and its power spectrum (right panels): around (a) x ~ 2300A C and (b) x ~ 2350A C . The 
portions where the power is strong are consistent with the linear theory of (a) the Buneman instability and (b) the ion-acoustic instability, 
respectively. 

Thus, it can be concluded that the filamentary structure in the foot region is generated by the ion beam-Weibel 
instability. The strong magnetic field generated by the instability would contribute to the thermalization of the 
incoming ions immediately upstream of the ramp (see Fig.[7])- 

3.3. Downstream temperature and jump condition 

In the downstream region, we obtain a temperature ratio of T e /T; ~ 0.38 from Fig. [7j Thus, the ratio is significantly 
smaller than unity, although it is still much larger than those observed in several SNRs [e.g., T c /T p < 0.07 in SN1006; 
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Fig. 12. — (a) Ion number density, (b) magnetic field B z , and (c) current density in the x-direction, J x , around x ~ 2375A C . (d) Power 
spectrum of J x calculated in the region 2362.5 < x/\ B < 2385.5. 

iGhavamian et al. I (|2002t ).] 

Figure [14] shows the kinetic energy distributions of the electrons and the ions in a rectangular downstream region 
(2624 A fi < x < 2656 A e ). Both d istributions are fitted very well with the (3D and relativistic) Maxwellian distributions 
(e.g., lLandau fc Lifshitz lfl980h 

/( 7 )d7 cx 7( 7 2 - 1) 1/2 exp(-- T mc 2 /k B T)dj (19) 

with temperatures k B T c jm c c 2 = 0.14 (for the electrons) and k-QTi/m c c 2 — 0.42 (for the ions) for E\a n /m e c 2 < 2. 
These temperatures again give a low temperature ratio of T a jT\ ~ 0.33. When the upstream bulk kinetic energy 
is completely dissipated into thermal energy and the electron and the ions are in equipartition, the temperature is 
given by k B T jm^c 2 ~ 0.32 for both species. On the other hand, when the electrons and the ions are thermalized 
separately, their temperatures are kBT e /m e c 2 — 0.021 and k B T\lm c c 2 = 0.63, respectively. The ion distribution has 
a suprathermal tail for Eki n /m c c 2 > 2. As shown in the next subsection, these suprathermal ions originate from 
the reflected ions in the foot region. In contrast, neither a suprathermal tail nor an accelerated population is clearly 
observed in the electron distribution. 

Using the above results, the shock jump conditions are calculated as follows. In the shock rest frame, the upstream 
flow velocity V\ and the downstream flow velocity V% are given by V\ = 0.33c and V% = 0.084c, respectively. Thus, we 
have Vi/V 2 ~ 3.9, A^ 2 /Ai - 4.1, and (k B (T c + T^/rrii) 1 / 2 - 0.14c, where N x {= n c0 ) and N 2 arc the number densities 
in the upstream and downstream regions, respectively. On the other hand, the MHD Rankine-Hugoniot relations 
(e.g.. iTidman fc Krall lH97ll ) give V1/V2 = N 2 /N 1 - 4 and (k B (T c + T i )/m i ) 1 / 2 - 0.15c in the high Mach number 
limit. Hence, the simulation results agree very well with the MHD Rankin-Hugoniot relations. Even although, in the 
downstream region, the magnetic field reaches ~ 15 times the upstream fields [see Fig. EJb)], the plasma beta is still 
high ((3 ~ 25) and the magnetic pressure is negligible for the jump condition. 

3.4. Acceleration of reflected ions 

As Fig. [T4l shows, a fraction of the ions are slightly accelerated to E^ n /m e c 2 ^3 — 6 (measured in the downstream 
frame) and form suprathermal populations. Figure [TST a) shows the trajectories of two typical accelerated ions (red and 
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Fig. 13. — (a) The linear growth rate and (b) wavenumber of the most unstable mode of the electromagnetic instability with the wave 
vectors parallel to the background magnetic field (i.e., in the ^/-direction) (solid curves). The real frequency is shown by the dotted curve 
in the panel (a), which shows that the mode is a purely growing mode. Those for the ion beam-Weibel instability are shown by the dashed 
curves (for V BX ) and dot-dashed curves (for V BZ ), respectively. 



blue curves) together with that of a non- accelerated ion (green curve). It is clear that accelerated ions are reflected 
at the shock front (i.e., the ramp) and go around the upstream region, whereas non-accelerated ions are directly 
transmitted downstream. Figure [TBTb) shows the kinetic energy history of the three ions measured in the downstream 
frame. It reveals that the kinetic energies of the ions increase while they are in the upstream region. This is simply 
acceleration by the motional electric field E z = —V X B V jc in the downstream frame while they are gyrating in the 
foot region (lAuer et al. lll97lHPhillips fc Robson Ill972t ) after specular reflection at the shock front (|Paschmann et al~l 
1982; Goslin g et al. 11198$ iSchwartz et al. 111983ft . As Fig. Wis) shows, the ions are accelerated when propagating in 
the — z-direction. This process can be understand more clearly in the upstream frame where there is (essentially) only 
a background magnetic field and no motional electric field. Figure I15f d) shows the same kinetic energy histories as 
Fig. 115( b) but measured in the upstream frame. It shows that the ions gain energy at the reflection and subsequently 
their kinetic energy remains almost constant. Thus, the ion acceleration is simply due to reflection at the ramp. 

3.5. Currents and magnetic field 

Figure [16] shows each component of the current density and the magnetic field. The upstream background field, 
which is in the y direction, is compressed in the shock transition region as in the ID simulations, although it fluctuates 
considerably in the present case. As discussed above, many current filaments exist in the foot region of the shock in J x 
and J z ] the filamentary structures observed in the ion number density in Fig. [2] indicate the presence of these current 
filaments. The filaments generate a magnetic field in the same way as a Weibel-mediated shock in unmagnetized 
plasmas except that in the present cases B x and B y components are generated by the current filaments in addition to 
the B z component because the background field deflects the particles in the z-direction and the current filaments can 
have a J z component as well as a J x component. 

Figure [171 shows the magnetic field strength normalized by the upstream background magnetic field, \B\/Bq. There 
are some strong, highly tangled magnetic fields in the transition region (\B\/Bq ~ 40) and the downstream region 
(\B\/Bq ~ 15, which is much larger than the magnetic field strength when merely compressed ~ 4). Figure IT51 shows 
enlargements of the current density in the z direction, J z , and the magnetic field strength \B\ in a rectangular area 
in the downstream region (2550A C < x < 2700A C and 50A o < y < 200A o ). It is evident that the downstream tangled 
magnetic field is mainly generated by the current filaments in J z . The filaments have typical sizes of ~ 2 — 4Ai, where 
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Fig. 14. — Kinetic energy distributions of the electrons and the ions measured in the downstream frame (solid histograms) normalized 
by the electron rest mass energy in the downstream region within 2624A C < x < 2656A C . The dashed and dotted curves are the Maxwellian 
distributions with temperatures of k^T e /m c c 2 = 0.14 for electrons and k^T\/rn c c 2 = 0.42 for ions, respectively. 

A; = (mi/mo) 1 / 2 A e is the ion inertial length, which is slightly larger than the filament size in the foot region. This 
can be explained by current filaments coalescing downstream of the foot region. Some of the current filaments have a 
complex coaxial structure in that they are surrounded by return and anti-return currents. 

Coalescence of current filaments that carry a current in the cc-direction in the foot region is inhibit ed due to the 
dimensionality of the simulation, while current filaments in the z-direction can merge with each other fc.f. iMorse 1119711 ; 
iLee fc Lampe |[T973l : lKato~l l2005). Furthermore, the current filaments in the z-direction are not affected by instabilities 
in the current direction, such as the kink instability. Therefore, the current structures in the foot and downstream 
regions may still differ in three dimensions. 

3.6. Early evolution 

It is interesting to note the evolution of the system before the effect of the background magnetic field becomes 
significant (i.e., t < Q,^ ~ 2000W",, 1 ). In this period, the plasma is effectively unmagnetized and another kind of 
shock, namely an unmagnetized shock, appears. Figures [T9l and BUI respectively show the evolution of the ions in the 
x-u x phase space and that of the ion number density. At a very early time (co pe t = 500), the incoming ions and 
the ions reflected by the wall at x = 3000A G form a counterstreaming beam system and the number density in the 
overlapping region simply becomes twice the upstream number density. The two populations then start to interact 
(1000 < uj pc t < 2000) and the density on the right side of the overlapping region increases. As shown below, this 
interaction is due to the magnetic field generated by the ion beam-Weibel instability in the overlapping region that 
deflects the ions; this provides a kind of dissipation mechanism. At a later time (3000 < ui pc t < 4000), the effect of 
the background magnetic field becomes important and incoming ions start to accumulate around x ~ 2775A C . Most 
of the ions that are initially reflected by the wall gyrate back downstream due to the background magnetic field by 
uj pc t — 3000. Instead, at oj pc t — 4000, another reflected ion population appears around x = 2700A o ; the ions in this 
population have been newly reflected at the 'ramp' (x ~ 2750A C ). Thus, the structure changes from an unmagnetized 
shock to a magnetized one around that time. 

Figure [21] shows the ion number density, the x-component of the current density J x , and the z-component of the 
magnetic field B z at uj pc t = 2000. It is evident that many current filaments are generated in the interacting region and 
generate a strong magnetic field around themselves. A filamentary structure is also observed in the number density 
corresponding to the current filaments. These current filaments are generated by the ion beam-Weibel instability 
between the counterstreaming ion populations. The generated magnetic field provides the dissipation mechanism for the 
ions to form the unmagnetized shock. Figure [22] shows u x -u y plots in the following three regions: 2700 < x/X c < 2800, 




Fig. 15. — Histories of two accelerated ions (red and blue curves) and one non-accelerated ion (green curve): (a) Trajectories on the 
evolution of E z . (b) Kinetic energies of the ions measured in the downstream frame, (c) ^-coordinates of the ions, (d) Kinetic energies of 
the ions measured in the upstream frame. 
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Fig. 16. — Left panel: each component of the current density (a) J x , (b) J y , and (c) J z normalized by n e oec. Right panel: each component 
of the magnetic field (d) B x , (e) B y , and (f) B z normalized by Bo. There are numerous current filaments in the shock transition region 
and the downstream region and they generate the magnetic field. 
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Fig. 17. — Magnetic field strength normalized by the upstream background magnetic field, \B\/Bq. There is a strong, highly tangled 
magnetic field in both the shock transition region (\B\/Bq ~ 40) and the downstream region (\B\/Bo ~ 15). 
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Fig. 18. — (a) Current density J 2 normalized by n e oec and (b) magnetic field strength |B| normalized by the upstream field Bo in the 
downstream region 2550A C < x < 2700A C and 50A C < y < 200A C . There are numerous current filaments carrying currents in the z-direction 
and they mainly generate the downstream tangled magnetic field. Some of the current filaments have a coaxial structure. 
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Fig. 19. — x-u x phase space plots of the ions in the early evolution. Phase space plots at (from top to bottom) Lo pe t = 500, 1000, 1500, 
2000, 3000, and 4000. 

2800 < x/\ c < 2900, and 2900 < x/X c < 3000. The ions are almost completely isotropized and form a ring-like 
distribution in the u x -u y plane far downstream (2900 < x/X c < 3000). Since the beam-Weibel instability generates 
only the z-component of the magnetic field in this 2D configuration, the magnetic field deflects the ions only in the 
u x -u y plane. However, in three dimensions, the generated magnetic field can have the y-component and the ions can 
be deflected in all directions, resulting in three-dimensional dissipation. 

The dissipated ions form the downstream region of the shock and the shock structure propagates upstream at an 
almost constant speed, as shown in Fig. Q] for 1500 < u> pc t < 2000. In the present simulation, this Weibel-mediated 
shock disappears when the effect of the background magnetic field becom es significant at later t imes. However, this 
shock can propagate steadily upstream when there is no background field (|Kato fc Takabe Il200l . 



4. DISCUSSION 
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Fig. 20.— Time evolution of the ion number density for uj pe t = 500, 1000, 1500, 2000, 3000, and 4000. 



In the previous section it was found that the filamentary structures in our simulation are generated by the ion beam- 
Weibel instability in the foot region. Similar structures have been found in the foot or overshoot region in 2D PIC 
or hybrid simulations with lower Mach numbers. There are two different causes for these structure s: the emission of 
whistler waves at the ramp (jKrauss-Varban et al. 1119951 IHellinger et al. |[2007t ILembege et al. 1120091) and the emission 
of Alfven waves due to the A lfven ion cyclotron instability, resulting in a structure called rippling (|Winske fc Quest 1 
[19881: ILowe fc Burgessl feOOS) . Since these are associated with waves generated at the ramp or the overshoot and which 
then propagate upstream, these processes could als o be related to t he filamentary structures in our simulation. The 
wavenumber range of the whistler wave is given by (|Ichimarulll973l ): 



4m c n,,< 



kc 

pc 



«1 



where n c is the local electron number density at which the whistler wave exists and n c o and uj 
far upstream. The left-hand side is ~ 0.133 in our simulation and if we take n e /n o ~ 2 [Fig 

[Fig. IT2Td)] as typical values, the above condition is (marginally) satisfied. However, if the structure is related with 
(standing) whistler waves, its group velocity must be greater than the shock speed (in the upstream frame): 



(20) 

are defined for the 
:)] and kc/uj pc ~ 1 
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where B$ and uj cc are defined for the far upstream, whereas B is the l ocal value. [Here, we neglec t the dependence 
of the propagation angle and hence it is just a necessary condition; c.f., iKrauss-Varban et al. I (|1995T ) .] This condition 
can be rewritten as 



M A < 2 



-Bo \ n e0 



kc 



(22) 



J y>c 



In our simulation, the right-hand side is ~ 10 — 30, whereas the left-hand side is ~ 130. Thus, the filamentary 
structure observe d in our simulati o n doe s not originate from whistler wave emission. On the other hand, the result of 
the simulation bv ILembege et al"~l (|2009( ) with m\/m c = 400 satisfies this condition if B/Bo ~ n c /n c o; the right-hand 
side is ~ 2fcc/w p i ~ 14.6, whereas the left-hand side is ~ 4.93. Assuming kc/uj pc ~ 1 and B cx n Cl a rough condition 
for the Alfven Mach number is given by 

M A < 2(TOi/m ) 1/2 . (23) 

This gives Ma < 11 for the mass ratio in our simu l ation (m-Jm P , = 30) a nd Ma < 86 for the real mass ratio of 
mi/m e = 1836. The rippling (Wins kefc Quest~lll988ULowe Sz Burgess 1120031 ) cannot be the cause of the filamentary 
structure in our simulation because (1) the rippling develops in the overshoot region, whereas the filamentary structure 
in our simulation develops in the foot region; (2) the wavelength observed in our simulation (A ~ 8A C ) is smaller than 
that of Alfven waves (Xa > 3A; ~ 16A e for mi/m e = 30); and (3) the Alfven speed in the f oot is smaller than the shock 
speed. Nevertheless, comparison of Fig. [16] with the figures in iWinske fc Quest"! (Qj)88) reveals that the structures 
of the magnetic field and the number density at and behind the overshoot are similar to each other. Therefore, the 
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21. — (a) The ion number density, (b) the ^-component of the current density, and (c) the z-component of the magnetic field at 
u pC D — 2000. It is evident that many current filaments exist and generate a strong magnetic field, which eventually isotropizes the incoming 
ions. 
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rippling mechanism may also operate there in our simulation. iLembcgc et al. (2009]) showed that these structures 
can also be affected by the dimensionality of the simulations in 2D simulations; the structures when the background 
magnetic field lies in the simulation plane may differ from those when it is perpendicular to the p lane. 

Electron acceleration was not observed in our simulation, whereas lAmano fc Hoshinol (|2009ft observed a kind of 
electron acceleration in a perpendicular shock in their 2D PIC simulation, which used similar parameters to ours. The 
main reason for this is considered to be the different direct ions of the upstream back ground field: in our simulation 
it lay in the simulation plane, whereas in the simulation by lAmano & Hoshino (2009|) it was out of the plane. This 
again demonstrates that dimensional ity can affect the results, even for 2D simulations. 

The ion-ion streaming instability ( Stringer! Il964t lOhnuma fc Hattalll965t iForslund fc Shoiikl 11970) . which has a 
wavevector that is highly oblique to the streaming direction, can be driven by the interaction between t he incoming 
ions and the reflected ions streaming upstream in the foot region and it ca n contribute to ion heating (|Auer et al. I 
119711 : iPapadopoulos et al. 1119711 : IWu et alTI 11984 lOhira fc Takaharal 12008). [This instability can also be driven in 
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Fig. 22.— u x - u y plots of the ions at LU pE t = 2000 in three regions: (a) 2700 < x/\ c < 2800, (b) 2800 < x/X c 
2900 < x/\ c < 3000. In (c), the ions are almost isotropized in the u x - u y plane, forming a ring-like distribution. 



< 2900, and (c) 



front of electrostatic shocks in two dimensions (jKato fc Takabe 1I2010D .] However, this instability was not observed in 
our simulation (see the right panels in Fig. IPTl). This is due to the high temperature of the reflected ions streaming 
upstream, Ir_, in our simulation. According to lOhira~fc Takahara (2008]), the ion-ion streaming instability is efficient 
for wavenumbers k > kr> e , whereas it is damped for k > kui due to the thermal motion of the ions. Therefore, ku c < kui 
is a necessary condition for efficient growth of the instability. In the present case, the two ion populations have different 
temperatures (Fig. [7]) and the Debye wavenumber of the reflected ions, /cd,r_ , should be used to evaluate the thermal 
damping effect. Thus, the condition for effective growth of the ion-ion streaming instability is fcoc < &d,r_ or 



Tr- UR- 



(24) 



The right-hand side of this equation is always smaller than unity, whereas the left-hand side is larger than unity in 
our simulation, as shown in Fig. [8l For example, at x = 2375A C , the ratios are Tr_/T c = 4.2 and riR-/n c = 0.32, 
respectively (see Table [J). Hence, the above condition is not satisfied and this would explain why the instability does 
not grow in our simulation. To examine the effect of Tr_ more quantitatively, we performed linear analysis with the 
parameters obtained from the simulation at uj pc t = 8000 and x = 2375 A c , but changing Tr_ in the same way as in 
Fig. IH1 Figure [53] shows the linear growth rates, the wavenumbers, and the angles between the wave vector and the 
streaming direction of the most unstable mode of the ion-ion streaming instability as functions of Tr,-. It shows 
that the instability depends strongly on Tr_. In particular, for Tr_/T > 0.35, it does not grow at all. Thus, the 
temperature of the reflected ions (Tr_) is important for growth of the instability in the foot region as well as the 
Buneman instability. 

Finally, we mention the possibility of generating magnetized shocks in experiments . Pres ent large-scale laser facilities 
can generate collisionless plasma flows at speeds of ~ 1000 km s _1 (iTakabe et al. I [200 81). Thus, if magnetized colli- 
sionless plasmas flowing at this velocity can be generated in laboratories, it should be possible to perform experiments 
on magnetized shocks. Table [2] shows the required background magnetic field strengths Bq for several sigma values 
and for a number density of n e o = 10 20 cm~ 3 together with the corresponding ion gyration time, T g , and ion gyro 
radius, r g . In this table, the gyration time and the gyro radius are calculated for an ion mass of 1836m e and a flow 
velocity of 1000 km s _1 . Sigmas of a = I0~ 3 or a = I0 -4 give achievable values for present large-scale laser facilities. 



TABLE 2 

Magnetic field strength required for experiments 



a 




M A 


Bq (G) 


T g (s) 




Tg (m) 




10" 


2 


10 


4.6 x 10 5 


1.4 x 10" 


ll 


2.3 X 10- 


4 


10" 


3 


32 


1.4 x 10 5 


4.7 x 10- 


-9 


7.3 x 10" 


4 


io- 


4 


100 


4.6 x 10 4 


1.4 x 10- 


-8 


2.3 x 10- 


3 


10" 


5 


320 


1.4 x 10 4 


4.7 x lu- 


-8 


7.3 x 10- 


3 



5. CONCLUSION 

We performed a 2D PIC simulation to investigate collisionless shocks propagating in weakly magnetized electron-ion 
plasmas at nonrelativistic speeds with a sigma of cr = 10~ 4 . We showed that current filaments are generated within 
the foot region by the ion beam-Weibel instability and that they generate magnetic fields in the same manner as for 
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Fig. 23. — (a) The linear growth rate, (b) wavenumber, and (c) angle of the most unstable mode of the ion-ion streaming instability 
as functions of the temperature of the reflected ions streaming upstream, Tr _ , with the other parameters obtained from the simulation at 
Wpet = 8000 and x = 2375 A e , as in Fig. [9] The arrow indicates the initial value of Tr_ in the simulation. It is clear that the instability 
significantly depends on Tr_. 



Weibel-mediated shocks in unmagnetized plasmas. The magnetic field strength generated by the current filaments is 
comparable with or even stronger than those of the compressed background magnetic field. Therefore, these filaments 
and their associated magnetic fields, which cannot be analyzed by ID simulations, are important in the formation 
of collisionless shocks in weakly magnetized plasmas. There are current filaments in the downstream region and 
they generate a tangled magnetic field that is typically 15 times stronger than the upstream background field. The 
thermal energies of electrons and ions in the downstream region are not in equipartition and their temperature ratio 
is given by T c jT\ ~ 0.33. We found a fraction of the ions were slightly accelerated on reflection at the shock, whereas 
significant electron acceleration was not observed in our simulation. The simulation results agree very well with the 
Rankine-Hugoniot relations. It was also shown that electrons and ions are heated in the foot region by the Buneman 
instability (for electrons) and the ion-acoustic instability (for both electrons and ions). However, the growth rate of 
the Buneman instability was significantly reduced from typical growth rates of this instability because of the relatively 
high temperature of the reflected ions. For the same reason, ion-ion streaming instability did not grow in the foot 
region. 
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